Linewidth of a polariton laser: Theoretical analysis of self-interaction effects 
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Polaritons in semiconductor microcavities can experience a Bose-Einstein condensation experi- 
mentally detectable in the spectrum of the emitted light. Scattering with noncondensed particles 
as well as self-interaction in the condensate provoke phase-diffusion limiting the coherence of the 
polariton condensate. We present a theoretical analysis of self-interaction effects on the lineshape 
of the emission from a polariton laser. Our calculations, for CdTe microcavities, show that there is 
an optimum pump at which the linewidth of the emitted light is reduced down to 1/j.eV. 
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A Bose-Einstein atomic condensate is a source of co- 
herent matter- waves (matter laser), in the same way as 
a photonic laser is a source of coherent light. Another 
system capable of undergoing a similar condensation is 
that of polaritons formed by the strong coupling between 
quantum well excitons and confined photons in a semi- 
conductor microcavity. Polaritons behave as composite 
bosons at densities below the saturation density, as con- 
firmed by recent experiments that include the observa- 
tion of stimulated scattering, and parametric amplifica- 
tion and oscillation [1,2]. A polariton Bose-Einstein con- 
densate is a matter laser that can be optically pumped, 
and experimentally detected by the emitted light [3-5]. 
The growth of semiconductor microcavities with new ma- 
terials, such as II- VI compounds, or GaN, opens great 
possibilities, due to the strong stability of the exciton 
in these systems [6,7]. In particular, a recent calcula- 
tion [8] shows that in high quality CdTe microcavities, 
huge occupation numbers can be achieved in the polari- 
ton ground state at densities well below the saturation 
density, i. e., at densities at which polaritons can be 
described as interacting bosons. Under these conditions, 
the system would be unstable to symmetry-breaking, and 
thus, to the formation of a BEC. 

In this work, we establish experimental signatures of 
a polariton laser. We consider a microcavity pumped by 
a continuous nonresonant laser, so that, the polariton- 
polariton scattering is fast enough to overcome the radia- 
tive losses and the system is able to relax to the ground 
state. At densities larger than a given threshold for BEC, 
the lowest energy level shows a macroscopic occupancy, 
and the system becomes a continuous polariton laser. 
The standard theory of photon lasers [9], i. e. non- 
interacting systems, predicts a very narrow linewidth, 
inversely proportional to the number of condensed par- 
ticles, T NI oc 1/Nq. However, polaritons interact with 
each other through some potential V and new physics is 
involved. Self-interaction in the condensate provokes a 
process of phase-diffusion that is determined by the en- 
ergy scale VNq [10-12]. When VNq is comparable to 
KT NI , this process would increase the linewidth as the 



number of condensed bosons increases, a behavior that 
is the opposite to that of a photon laser. We introduce 
here a self-consistent framework to include the two ef- 
fects described above. We show that there is an optimum 
pumping range to get an extremely narrow linewidth of 
the polariton emission. For CdTe microcavities, we find 
that the linewidth can be reduced down to lfieV. 
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FIG. 1. Phase space for the scattering of polaritons in a 
semiconductor microcavity. Upper and lower dashed lines 
represent the bare photon and bare exciton, respectively. The 
continuous line is the polariton dispersion that results from 
the strong exciton-photon coupling. 



Fig.l depicts schematically the lower energy branch of 
the spectrum of a semiconductor microcavity for the case 
of zero detuning [14]. Below the minimum bare exciton 
energy, eg, excitons and photons are strongly coupled. 
The ground state (k = 0) lies at the energy ftp/2 below 
eg, where Qp is the polariton splitting. At energies above 
eg, the polariton branch merges abruptly into the bare 
exciton dispersion, and the exciton-photon coupling can 
be neglected. We consider the phase space as divided in 
two regions: the lower energy states, labelled as lower 
polaritons (LP), and the states above eg, treated as bare 
excitons. The small density of states of the LP, when 
compared to the exciton one (plp ~ 10~ 4 Px), drastically 
reduces the threshold density for the formation of a BEC 
in the LP part of the spectrum. This allows us to simplify 
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the problem of the polariton dynamics by considering the 
exciton states as a thermalized reservoir [8] . On the other 
hand, microcavity polaritons are 2-D quasiparticles that 
can undergo a Kosterlitz-Thouless transition in the ther- 
modynamic limt [13]. Finite size effects leads to a local 
transition to BEC, so that we quantize the LP levels ac- 
cording to a scale determined by area S of the spot of 
the pumping laser. 

In this work we are not analyzing the true quantum 
state of the polariton condensate, but we are interested 
in the spectrum of the emitted light, that can be calcu- 
lated from the adequate correlation function. The steps 
of the theoretical analysis are: first, to obtain an equation 
of motion for the density matrix, and later to compute 
correlation functions and emitted intensities. 

Equation of motion for the density matrix. The main 
mechanism for the relaxation to the lower energy states 
is the scattering of two excitons [15], in which one of 
the final states is a LP. We label this scattering pro- 
cess exciton-polariton (XP) scattering. This relaxation 
process creates a non-equilibrium polariton distribution 
that evolves towards a Bose-Einstein condensate as ex- 
citon density is increased [8]. Our Hamiltonian is the 
sum of three terms describing the bare-exciton and LP 
dispersions (Ho), the XP scattering (Hxp), and the self- 
interaction in the condensate mode (Hsi). 

Ho = Y, e k Pa W +Y, e k b l bk ' 

k k 

H X p = y^<4-Ffc + h.c, H S i = Valala a a . (1) 

k 

The index k is quantized according to S. a\ 
(e k p ), bl (ej?) arc the creation operators (energies) 
of the LP, and bare excitons, respectively. F% — 
Y,k 2 ,k 3 M Vk ^' k 3M b k 2 b k3 b k 4 describes the scattered ex- 
citons, with VkxMMM being the polariton-polariton in- 
teraction in the kap « 1 limit (see [16] for an ex- 
plicit expression). V = Vo,o,o,o is the self-interaction in 
the ground state. We have neglected polariton-polariton 
scattering that involves more than one lower polariton, 
and polariton-phonon scattering, because they are much 
slower than the process depicted in Fig. 1 , and only pro- 
duce energy shifts [8,16]. 

Since we are mainly interested in the evolution of LP, 
we trace out the reservoir degrees of freedom in the den- 
sity matrix operator x> and define the reduced density 
matrix s: 

s(t) = Tr R { X (t)}, (0 LP (t)) = Tr LP {0 LP (t) S (t)}, (2) 

where Trp, Tr^p represent the trace over the exciton 
reservoir and the LP, respectively, and Olp is any func- 
tion of LP operators. 

We describe the exciton reservoir by a thermalized 
Maxwell-Boltzmann distribution. This approximation 



is justified by the fast exciton-exciton scattering within 
the exciton reservoir. Moreover, this assumption is sup- 
ported by a recent experiment [17] on the evolution of 
the polariton distribution when pump-power or temper- 
ature is varied. The total density matrix operator can be 
approximated by: 

x (t) » s(t) ® f R (t) = s(t) ® e E»("-- £ !) i l 6 */ fc B T - ) (3) 

where /j, x , T x are the chemical potential and temperature 
in the exciton reservoir, respectively. 

The time evolution of the density-matrix is calculated 
in the interaction picture. Up to the lowest order in 
Hxp + Hsi, s evolves as: 

(a k al's'-al's'a k )(F k F^) + (a{a q 's'-a q 's'a{)(FlF q ) + 
(s'al'ak - a k s'al')(F^F k ) + (s'a q 'a\ - a\s' a q ')(F q ' Fl) , (4) 

where all operators are in the interaction picture with 
respect to Hxp + Hsi, and the primes indicate time de- 
pendence on t' . 

Eq. (4) becomes simpler when LP operators at time 
t are translated to the previous time t' by a\(t) = 
a\(t')e lf -k F '(t-t')/^ _ Moreover, we are interested on the 
steady-state regime, i. e., in the limit to — > — oo, so that 
a Markovian approximation is well justified. Under this 
approximation, we take the LP operators in Eq. (4) out 
of the time integration and obtain the following master 
equation: 

ds iV T t 2 2l V- ( W k n , t t s 

dt = Jt i s,a ° a ° \ ( ~K Sfl fe -«fc a fe s )+ 

k ^ 

k -(a k sa\ - a\a k s) + h.c.\ . (5) 

The rates W™ and W^ ut are easily evaluated with the 
thermalized exciton distribution [8]: 

w Mout )= J2 A\V k ,k2^ M \ 2 {l + Nl 2 )Nl 3 Nl 

k2,k 3 ,ki 

K4 P + 4 2 - 4 3 - 4 J, 

W k ut = £ 4|U fc , fc2 , fc3 , fc 4| 2 ^ 2 (H-iV^)(l + ^ 4 ) 

k 2 ,k 3 ,k 4 

s(4 p + 4 2 -4 3 -4J, (6) 

where are the occupancies in the exciton reservoir. 
The imaginary parts of the time integrations have been 
neglected because they lead to energy-shifts irrelevant 
for the process of BEC and phase-diffusion that wc con- 
sider here. In Eq.(5), we have also included a term (not 
appearing in Eq.(4)) [9] that accounts for the radiative 
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losses with a rate Tk = Cfc p /r p h, where C^ p is the pho- 
ton weight in the polariton wave function, and r p h is the 
lifetime of the photonic mode confined in the microcavity. 

Eq. (5) is the keynote of our analysis. In particular, it 
allows to calculate the evolution (which does not depend 
on the self-interaction) of the LP occupation numbers: 

j t {a\a k ) = Wr({a\a k ) + 1) - (W^ + r fe )<4^>- (7) 

We also describe self-consistently the evolution of the pa- 
rameters n x , T x , that describe the exciton reservoir, by 
deriving the corresponding rate equations as described 
in [8]. From the steady n x , T x , we get the rates Wq™, 
Wo"*, that correspond to the scattering to the ground 
state. After that, we have found all the parameters in 
the master equation (5). 

Emission spectrum. Once we have a self-consistent de- 
scription for the time evolution of the density matrix, we 
can undertake the task of computing expectation values 
of magnitudes experimentally measurable. In particular 
the emission spectrum can be obtained from the two-time 
correlation function: 

1(e) = -Re / (al{T)a {Q))e- lf - T l n dT. (8) 
t Jo 

Without self-interaction in Eq. (5), the application of 
the quantum regression theorem would lead trivially to 
Lorentzian line shape of the spectrum: 

(aUr)a om ^ = N e^ N ^, T»' = ^±?±. (9) 

Decoherence has the usual aspect of a single-particle 
noise corrected, in the denominator, by the number of 
particles in the condensate. However, inclusion of the 
self-interaction term changes this result dramatically. We 
use a well-known method in quantum optics that allows 
to include exactly the effect of the self- interaction [9]. 
The two-time average is expanded as a sum over a set of 
auxiliary functions C n (r) : 

(al(r)a (0))=J2V^C n (T), 

n 

C n {T) = {e tHoT l h \nXn-l\e- lHoT l h a Q {Q)), (10) 

with \n > being the number representation of the k = 
mode. One can use the quantum regression theorem, to 
show that the functions C„ satisfy the differential equa- 
tion [10]: 

^Cn(r) = ^ p C n (r) + Ln, m C m (r), (11) 

m 

where the L„_ m are the coefficients governing the evolu- 
tion of off-diagonal one-time matrix elements: 



— < n- 1\s(t)\u > = ^£„, m < m- l\s(T)\m > . (12) 

m 

Using Eqs.(5) and (12), Eq. (11) becomes: 
^C n {r) = (-w+ (2n + 1) - w (2n - 1) + 

l -V{n -I- N )) C„(t) + w+ 2 v /n(n-l)C„-i(r) 

+ Wq ~ 2v/n(n + l)C„+i(T), (13) 

with w+ = W in /2 and w = (Wg ut + T )/2. In Eq. (13) 
we have taken the origin of energies at e$ p + VN , in or- 
der to compare the linewidths of the emission spectrum 
at different densities. Initial condition for Eq. (13) is 
obtained from C n for r = 0: 




where < n\s\n > is easily deduced from having a time- 
derivative equal to zero in the steady state. Instead of 
solving numerically the enormous set of Eqs. (13), the 
problem can be simplified, in the case No >> 1, by re- 
placing the index n by a continuous variable [18], obtain- 
ing a partial differential equation: 

—i-^ =„(«,„ + w +)—C(n,r) + 

d 

(2n(w -w+)+ w - + w +) —C{n, r) + (15) 
(W - w+) ^±!^ + _ No _ i)) C (n, r). 

In the limit HT NI >> VN n , (15) gives the adequate 
Lorentzian shape of the spectrum with a linewidth T NI . 
In the opposite limit, VN n >> hT NI , an analytic solu- 
tion exists: 

C(n,r) w 4-V^e-^ +i % (n - No)T , (16) 
Wo 

predicting an asymmetrical lineshape: 

J( C ) „ ^ e + VN o) e -(e + VN )/(VN ) d(e + y No y 

V V No 

Results for a CdTe microcavity. We have solved nu- 
merically Eq. (15) for the case of a CdTe microcavity 
with Vlp = lOmeV, and zero detuning between the bare 
exciton and the photonic mode. The pump is assumed 
to add excitons at a a given rate, p x , and at a lattice 
temperature, 7l = 10K. This implies a very fast re- 
laxation by the exciton-phonon scattering, which is not 
always the case in experiments. The steady value of T x , 
however can reach 30 K, as explained in [8]. Lifetimes of 
the photon and the bare excitons are taken r p h = lps 
and t x — lOOps. The steady-state polariton density 
considered in our calculation is always below 0.3 times 
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the saturation density, which can be estimated to be 
6.7 x 10 n cm~ 2 in a CdTe microcavity [8]. The quan- 
tization length is 50fim, of the order of typica 1 excita- 
tion spot diameters. Fig. 2 gives No as a function of 
the pump-power. It shows a threshold for BEC (N > 1) 
around p x w 810 8 cm~ 2 ps~ 1 . 




FIG. 2. Full Width at Half Maximum of the emission spec- 
tra (continuous line) and ground state occupancy (dashed 
line) as a function of the pump-power p x (in 10 cm~ ps~ ). 
The optimum pump corresponds to p x « 10 9 cm~ 2 ps~ 1 . 



larger densities the linewidth of the polariton laser in- 
creases linearly with the number of condensed particles, 
until it reaches values comparable to the non-condensed 
emission. 

In conclusion, we have presented a theoretical method 
for the self-consistent calculation of the emission 
linewidth of a polariton BEC. Self- interaction sets an im- 
portant restriction for the coherence that can be achieved 
in this system. The dynamics described by our Eq. (5) 
would drive the system to a BEC provided a small coher- 
ent seed, i.e. a symmetry-breaking term, is included in 
the initial condition. Our conclusions are also relevant for 
the case of recent proposals in which the scattering mech- 
anism for the relaxation of polaritons towards the ground 
state is different than the XP scattering considered here, 
as for instance the case of the proposed electron-polariton 
scattering in doped microcavities [19]. 

We thank H. Kohler and J. Fernandez- Rossier for fruit- 
ful discussions. Work supported in part by MCYT of 
Spain under contract MAT2002-00139 and CAM under 
contract 07N/0042/2002. 
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Right: Continuous increase of the linewidth shown in the 
spectra for p x = 1.5(inner), 2.0, 3.0(outer) in units of 
W 9 cm- 2 ps- 1 . 

For densities below the condition N n V ~ HT NI , 
the calculated lineshape is Lorentzian with minimum 
linewidth of the order of 1/ieV. After condition NqV w 
KT NI is reached, there is a transition from the laser-like 
decoherence to the self-interaction broadening as shown 
in Fig. 3 (left). For larger pump-powers, the emission 
is asymmetric and is continuously broadened. Above 
p x = 1.110 9 cm _2 ps _1 , the numerical results are iden- 
tical to the approximation given in Eq. (17). In the 
evolution of the spectrum linewidth as a function of the 
pump-power, shown in Fig. 2, one observes that the op- 
timum pump is very well defined by the abrupt dip in 
the emission linewidth. This dip means an increase of 
the coherence of almost three orders of magnitude. For 
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